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Abstract 



We investigate the possibility of detecting CP-odd angular correlations in 
the various decay modes of the neutral Higgs boson including the modes of 
a ZZ pair, a W + W~ pair, or a heavy quark pair. It is a natural way to 
probe the CP character of the Higgs boson once it is identified. Final state 
interactions (i.e. the absorptive decay amplitude) is not required in such 
correlations. As an illustrative example we take the fundamental source of 
the CP nonconservation to be in the Yukawa couplings of the Higgs boson to 
the heavy fermions. A similar correlation in the process e + e~ — > l + l~H is also 
proposed. Our analysis of these correlations will be useful for experiments in 
future colliders such as LEP II, SSC, LHC or NLC. 
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I. INTRODUCTION 



The Higgs boson sector remains the most mysterious part of theories of electroweak 
unification. To probe this illusive sector one is well advised to keep an open mind. Most 
analyses of the Higgs sector, in the Standard Model or beyond, ignore potential CP violation. 
However, in many these models, even the lightest neutral Higgs boson can have interesting 
CP violating phenomena. In many models of CP violation, including the Standard Model, 
CP violation is a consequence of simultaneous existence of many coupling constants. Such 
requirement results in the presence of many coupling constants and/or loop suppression 
factors in the observables. The neutral Higgs sector is unique in the sense that a single 
Higgs boson coupling to a massive fermion is enough to manifest CP violation as long as the 
Yukawa coupling contains both scalar and pseudoscalar components. Therefore, in many 
ways the CP violating aspect may be the most interesting part of the Higgs boson physics 
beyond the the Standard Model once a neutral Higgs is identified. 

In a previous paper, we investigated the signatures of various CP-odd asymmetries 
in different polarized decay modes the neutral Higgs boson Among the interest- 

ing modes are the CP asymmetries in the event rate differences N(QlQl) — N(QrQr), 
N{W£W£) - N{W^Wr) or N{Z L Z L ) - N{Z R Z R ). As expected, the CP asymmetries will 
manifest themselves only after the final state interactions are taken into account. For the 
heavy quark pair, QQ, mode or the W + W~ pair mode, the relative energy of the final 
state charged leptons can be used as the polarization analyzers of the heavy quarks or the 
W's However for the case of the ZZ mode, the Z coupling to charged leptons is 

mainly axial-vectorial. Therefore the efficiency of using its leptonic energy spectrum as the 
polarization analyzer is suppressed by a factor of cv/ca = — 1 + 4 sin 2 9w ~ —0.08. 

In this paper we consider a different kind of signal of CP violation: the CP-odd correla- 
tions of final state momenta. Note that this is different from the more conventional T-odd 
signals which can be imitated by CP conserving final state interactions. No final state in- 
teraction is required to obtain the signals. Therefore, unlike the case of the CP asymmetry 
considered in Ref. 0, it is not necessary to cross a heavy fermion threshold to obtain a 
large signal in the W + W~ or the ZZ channel. For earlier discussions of CP-odd effects in 
different processes, see Refs. @-||. 

II. H -> ZZ 

Let us start with the ZZ mode whose CP-odd asymmetry does not translate into large 
lepton energy asymmetry. After the Z's decay there are three modes of final state: (1) 
/~/ + /'~Z' + , (2) l~l + qq, (3) qqq'q'. The purely leptonic modes are of course the easy ones in 
which to get a CP-odd correlation provided that we collect enough events. In the bigger 
samples of the semileptonic modes and the hadronic ones, it is difficult to tell the charge of 
the leading quark when it hadronizes into a jet. However, the most intriguing part of our 
result is that it is not necessary to identify charge to detect the CP-odd signal, an argument 
we will elaborate on later. This observation will be even more crucial in trying to decode 
CP-odd signals in the W + W~ or ti modes of the Higgs decay. Another important point is 
that in modes (1) and (3), the lepton pairs or the quark pairs do not have to be distinct. As 
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we will show later, the identical particle effect does not wash out the main CP-odd signal. 
We shall discuss these modes in order. 

Consider H — > ZZ — > If H is scalar, the two Z polarizations will be parallel; 

if H is pseudoscalar, the two Z polarizations will be perpendicular. The correlation between 
the polarizations of the two Z bosons translates into the correlation between the two planes 
of the l + l~ pairs. In fact this observation was used a long time ago to measure the parity 
of the 7T° through its decay into two photons ||. When the Higgs coupling is neither scalar 
nor pseudoscalar, one has a source of CP violation. We shall explore this case. 

The effective interaction for the ZZH vertex can be written as 

Chzz = \g H H[BZ v Z v + \DM H 2 e, vpa Z^Z^} . (1) 

where gu = {giMzl cos 6^), Z^ v is the field strength of Z boson and B and D are dimen- 
sionless form factors. For the Higgs boson in the Standard Model, D — 0, B — 1. The D 
term is CP-odd while the B term is CP-even. Simultaneous presence of D and B is CP 
violating. We have ignored higher dimensional CP-even operators in the above Lagrangian 
and a possible CP violating HZ'-f vertex which could also contribute to four fermions final 
states. Their effect should be negligible at the pole of the final state Z boson. In momentum 
space, these effective vertices give rise to the covariant amplitude: 

iM H ^z{p, v ),z(P',V) = igu{Bi] ■ rj' + (D/M 2 H )e( V , //, P, P')] . (2) 

Here 77, rf are the polarizations and P, P' are the momenta of the two Z bosons. We define 
6(77,77' , P, P') = e llua ^r) fl rj' l/ P a P' 13 , with the convention e i23 = 1- The helicity amplitudes in 
H — > ZZ are represented as 

M + , + = g H (B - iD0), M_,_ = g H (B + iD/3), M , = -^sl±|! , (3) 

in the H rest frame, with (3 2 = 1 — AM^/M H . Here the subscripts (+, — , 0) denote respec- 
tively the helicities (R,L, ||) of the Z bosons. Conservation of angular momentum implies 
only the above three helicity configurations for this decay process. Under CP transforma- 
tion, M + + — > M_ and M 00 M 00 . Thus, D is CP-odd, however Re D is CPT-even, 
while Im D is CPT-odd. The detail definition of CPT can be found in Ref . |§ . It is roughly 
the CPT symmetry without reversing the initial and the final states, that is, the symmetry 
reverses only the kinematic variables of the states according to their CPT property. The 
effect of Im D will give an asymmetry in the production of polarized states, 

N(Z L Z L ) - N(Z R Z R ) = (-/31m D/B)[N(Z L Z L ) + N(Z R Z R )} . (4) 

Its consequence in the energy asymmetry of the final leptons has been discussed in Ref. 

In this paper, we focus our attention only on the effect of the real part of D through 
the angular asymmetry of the final decay products. In the decay process H — > ZZ — > 
l~ (kJ)l + (k + )l'~ (p + ) , the angular correlation of the final 4 fermions is encoded in the 

3x3 density matrix 

Px = M\ \ My x , (no dummy summation), (5) 
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which is folded with the Z decay amplitudes to produce the event distribution. The CP-odd, 
CPT-even combinations are pi — p+ and p+ — p° + Po — Po while the ones that are CP- 
and CPT-odd are p\ — pZ and p° — p°_ — pZ + pg . Only the first two are relevant to our 
discussion below. 

Similar notation will be used for the modes of W + -W~. 

III. CP VIOLATING OBSERVABLES 

We will show observables which are related to the CP-odd correlations among the mo- 
menta, k + ,p_ and p + . To simplify our discussion, we start with the purely leptonic case 
and we arbitrarily label the lepton pair from one of the Z bosons with primes. The process 
Z — > 11 can be parametrized by the vertex 

i eu(l)jfj,(c L L + c R R)v(l) . 

Naively one may simply construct the CP-odd correlation 0( m = P- x p + ■ k_ in the H rest 
frame. It can also be written in a Lorentz invariant form 

O odd = -M H 1 e(p^,p + , fc_, k+) . 

However, as we will show later, the expectation value of this observable (p_ x p + • is 
proportional to cyCA where cy = \{cl + cr) and ca = \{— cl + cr) are the vector and axial 
vector couplings of the Z boson. Since the vector coupling of the Z boson to the charged 
leptons in the Standard Model is relatively small, this observable turns out to be rather 
unimportant. It can be understood as the consequence of an approximate symmetry when 
the vector coupling (or the axial vector coupling) is ignored completely, so that there is no 
distinction between l + and l~ as far as the Z boson is concerned. Therefore the differential 
decay rate is symmetric under two separate partial charge conjugation symmetries, C\ and 
Ci. Symmetry C\ interchanges l + and l~ while leaving Z /+ , / /_ unchanged; while the sym- 
metry C 2 interchanges / /+ and l'~ and leaving Z + , l~ unchanged instead. The usual charge 
conjugation operation C is the product of the the two C"s. It is easy to check that the 
correlation Q dd is odd under either C\ or C^. Therefore it has nonvanishing expectation 
value only when both cy and Ca couplings are nonzero. It is certainly more desirable to 
use an observable which is CP-odd but even under each C\ or C 2 - It is not too difficult to 
construct a quantity. A possibility is 

O'odd = (P- x P + ■ k-) [(p_ x p + ) ■ (k_ x %+)]. (6) 

It is easier to understand these angular correlations in the geometry of the reaction. A 
typical reaction is shown in Fig. 1. Under CP transformation, (E p _, Ek_,(f>) is transformed 
into (E p+ , Ek + , —(f)) . If one defines a plane using the cross product of one of the lepton pairs, 
say p- xp + , as the direction of the normal vector, the configurations with the momentum of 
the other lepton (fc_) coming out of the plane is the CP conjugate of the configurations in 
which k_ is going into the plane. Therefore the asymmetry, (sign(0 dd)) , can be interpreted 
as a kind of up-down asymmetry. More explicitly, define the polar and azimuthal angles for 
the two leptons in their respective Z boson rest frames to be 9, <fi for l(k-), and 9', 0' for 
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V{pJ). Here all the polar angles are defined relative to the same z axis, say the k_ + k + axis. 
The distribution depends on the relative azimuthal angle, therefore we can set azimuthal 
angle of I' to 0' = as shown in the Fig. 1. Then one can label the event configuration by 
the angles (cos#, cos#', 0). The energies E k _ and E p _ in the Higgs boson rest frame for the 
leptons I and V respectively are determined by their polar angles 9 and 9'. Under C\, 
parity and CP, the configuration transforms as follows, 



Cx 

c 2 
p 

CP 



(cos#, cos 9', i 
(cos 9, cos 9', ■ 
(cos 9, cos 9', . 
(cos#, cos 9', ' 



(C0S(7T — 9) , COS 9', <p + 71") 
(COS9, COs(7T — 9'), + 7r) 

(cos 9, cos 9', —(p) 
(cos(7r — 9), cos(tt — 9'),- 



(7) 



To understand observables such as (0' odd ), one can divide the azimuthal angle into four 
quadrants, I, II, III, IV. Since C\ transforms (I, II) into (III, IV), it is clear that the up-down 
asymmetry {sign(0 dd)) which corresponds to the angular integration of (I + II) — (III + 
IV) is odd under C\ or C 2 . On the other hand, an observable similar to (sign(0' odd )) can be 
constructed by the alternative angular integration (I+III)-(II+IV). 



IV. ANGULAR DEPENDENCE 

To calculate these asymmetries, we need the differential distribution: 

n(E p _,E k _,f 



dE p _ dE k _ dcf) 

A^EE Pxfx(8, 0)/-a(^ 0)fy(0, M-W, 0)44 • (8) 



hM A, A' 



We denote the helicity amplitude f\(9, 0) describing the process Z(\) —> l(h)l(—h) in the Z 
rest frame, where the spin projection of the Z boson along the z axis is A7i and the helicity of 
/ is specified by h = +(R) or —(L). One can relate the helicity amplitudes with the spin-1 
rotation matrix elements, f£ ~ d^ h e lX ^, 

fl(9,<j)) = (l±h cos 9)e ±l(f> /2 , f£ = hsm9/V2, (h = ±1) . (9) 

M is the normalization to be specified later. One can see from Eq. (HID) that in order to get 
the complex phase in angular distribution to expose the CP-odd, CPT-even effect, only the 
combinations pt — p+ and p Q + — p°_ + Pq — p$ contribute. That is consistent with our general 
discussion after Eq.(|5|). Note that in the H ZZ —> llll mode of identical final lepton 
pairs, the identical-particle symmetry implies n(E p _, E^^, 0) = n(Ek_, E p _, 0). Therefore it 
does not pose any restriction on the azimuthal, 0, angular distribution at all. The overall 
distribution can be divided into two parts, namely, the CP-even piece n and the CP-odd 
piece An, 

n(E p _,E h _,<f>) = n (E p _,E k _,(f)) + An(E p _,E k _ }( j)) (10) 
Using the form factors D, B and the couplings cl, c^, we obtain 
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An(E p _,E k ^,4>) = AfDM H 2 e(p„,p + , k-,k+) 

x h(ci + 4)(p_.fc_ +p+.A; + ) + 2c 2 L c 2 R (p-k + +p+.k-)] 
n (E p _,E k _,4>) = NB[(c\ + c R )p^.k- P+-k + + 2c\c 2 R p-.k + p + .kJ\ . 



(11) 



From now on, except for the Section IX, it is understood that we have dropped the Re prefix 
for the form factor D. We have only kept the linear piece in D, as a result of perturbation. 
Kinematically, we can express the e symbol and various scalar products in terms of angles, 



9,0,9', 



t(P-,P+,k-,k A 



\M 2 H M 2 z f3 sin 9' sin 9 sin. 



P- 


= 


p+ 


k + = 


P- 


•k+ = 


p+ 


■ k_ = 



M| 

16 

Ml 

16 

Ml 

16 

Ml 

16 



; 1 + p 2 ) (1 - cos 9 cos 9') + 2/3(cos 9 - cos 0')] - — £ sin 5 sin 0' cos . 

; 1 + /3 2 ) (1 - cos cos 0') - 2/3(cos - cos d')\ - —2- sin sin 0' cos . 

n M 2 

(1 + /3 2 ) (1 + cos 9 cos 0') - 2/?(cos + cos 9')\ + — ^ sin 5 sin 0' cos . 



1 + (3 2 ) (1 + cos 9 cos 0') + 2/3(cos 9 + cos 0')J + —2- sin sin 0' cos . 



(12) 



The covariant expression in Eq. (|TT| ) derived from the usual Dirac matrix calculation agrees 
with the result by the helicity amplitude method in Eqs.([|,|5],|],|[) with substitutions of 
Eq . (|12"1) . The differential CP-odd asymmetry then can be defined as 



•A u ,d. (B p _ , E k 



n(E p _, E k 



n(E p+ ,Ek + , 



n(E p _ , E k _ , 0) + n(E p+ , E k+ , 



An{E p _,E k 
no(E p _, E k _ 



(13) 



To simulate realistic detector acceptance, one has to use a Monte Carlo calculation based 
on the differential distribution in Eq.(|ll|). However, it is instructive to look at the overall 
asymmetry integrating over the full ranges of cos# and cos 9'. For the CP-odd numerator 
in A u .d. of Eq.([L3|), we have 



An(0) 



And cos 9d cos 9' 



MBDMfj i 
z (1 - 4^)2 



2-9 



z sin 20(c L + c R 



2 \2 



97T 2 

"64" 



1-2*) sm(f){c L -c R 



2 \2 



(14) 



with z = Mg/Ml. The first term in the bracket is due to the interference between M +j+ 
and M_ _ and is proportional to pi — p7. The second term in the bracket is due to the 
interference between M and M ±j± and is proportional to p\ — p°_ + Pq — Pq. For the 
CP-even denominator, we have 



no(0) = / no d cos 9d cos 9' 



MB 2 M% 



4-9 

9tt 2 



;i - 4z + 12z 2 + 2z z cos 20) (c2 + 4 



2 \2 



z(l-2z) cos0(c 2 -4) 2 



16 



(15) 
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The CP conserving part, no(0) has been calculated before []10 , |lT| . Our result agrees with 
Ref. Jllj]. A sign difference between ours and Ref. |HJ can be due to different definitions of 
6. The normalization M can be chosen to be 



Af = \%/[iiB 2 M A H {l 



4z + 12z 2 )(4 + 4) 2 ] 



(16) 



such that the distribution n(0) = n o (0) + An(0) is normalized to one after integration over 
0. Measuring these sin0 or sin 20 dependences in the event distribution will establish the 
CP non-conservation. To enhance statistics, we can look at the up-down asymmetry by 
integrating again over from to tt, 



A U a 



(So ~ SZ") n(<f>)d<f> _ D9irz(l-4z) 1 2(l-2z) 



fon(<f>)d<f> 



B 16(1 - Az + Viz 1 ) 



ci+ci 



(17) 



2(4 + 4 



gives zero 



After the integration, the sin 20 term of the combination c| + 4 
while the surviving sin0 term is relatively suppressed by the factor c| — c 2 R = —AcyCA as 
we anticipated before. A more realistic integrated asymmetry should be the one that can 
preserve the sin 20 term. That can be obtained by taking the difference in the integration 
over [0,7r/2] and the integration over [n/2,7i}. 



A,. 



( f^/ 2 f7T 

\J0 Jtt/ 



/2 



+ J? r/a -j£ /!1 )n(#ty 



Jo n(0)# 



D Az 2 {l-Az) 1 2 
'Bn(l-Az + 12z 2 ) 



(18) 



Fig. 2 shows A u .d.i A! u d versus Mh per unit of D/B. Following the above formalism, exper- 
iments can search for CP non-conservation, or at least, put a constraint on the parameter 
D. Note that, in the heavy Higgs limit z < 1, the sin0 contribution becomes relatively 
important because of the kinematic suppression factor z in the other contribution of sin 20 
in Eq.(Pj). 



V. GENERALIZATION 

We can generalize our formalism to the case that / and /' are different fermions of different 
couplings cl,cr and c' L ,c' R . The replacements in Eqs.( |IT| , |i4| -|T7j) are given by the following 
rules, 

4,4 2 '2 | 2 '2 

C L + C R ~ ¥ C L C L ' C R C R 5 
f) 2 2 2 '2 i 2 '2 

lc L c R ~^ C L C R ' C R C L ) 

(c| + 4) 2 -.(c| + 4)(4 + 4), 

(4 -4) 2 - (4 -4)(4 -4)- (19) 

This applies also to the case H —>■ ZZ — > l~l + qq and H — > ZZ — > qqq'q' . 

One of the main problems of using the hadronic or semi-hadronic decay modes of the 
Higgs boson in searching for CP violating signals is that experimentally it is impossible 
to identify charges of jets originating from hadronization of the partons. Therefore one 
essentially has to smear over the charge information. Naively one may think that it will 
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make CP information impossible to disentangle. This turns out not the case in the observable 
0' odd , which is even under both C\ and C 2 separately. Thus, in O' odd the contribution of an 
event and its Ci-conjugate event will add, while in O odd they will cancel. Therefore, charge 
identification is not necessary for the observable O' odd . 

To construct the observable, we revisit the definition of in the general case. Each 
pair of fermions defines a plane. Given two un-oriented planes one can define the angle 0o 
between planes to be a or 7r — a (0 < a < n by definition). These two choices are not resolved 
yet. However, if the common line of the two planes can be physically assigned a direction 
associated with one of the planes, then the two angles a and 7r — a can be distinguished. 
For example in the H — > ZZ decay, see Fig. 1, we can simply use the vector k_ + k + = kz 
to define a direction associated with the k-, k + plane. Then, using the right hand rule, one 
can rotate along the axis fc_ + k + and sweep the plane of p~,p+ toward the plane of k-, k + 
and define the resulting angle O £ [0, 7r]. It is important to notice that choosing p_ + p + to 

— * — # 

define the angle, instead of k_ + k + , gives the same result. (Therefore the two fermion pairs 
can be identical without smearing out the effect.) On the other hand, the azimuthal angle 
defined in Fig. 1 with full identification varies from to 2ir, and the relation between them 
is simply 0o =mod(0, 7r). So the event rate at 0o — o, (0 < a < % by definition) is the sum 
of the elementary rates n(0 = a) + n(0 = a + n), thus the dependences on sin0 in Eq.flUP 
and cos in Eq. ([15]) will be washed away. However, the sin 20 term survives to signal the 
CP violation. Of course, if one can identify the charges of the final fermions, as in the case 
of leptonic modes, then one can also decode the sin0 term as well. To detect the sin 20 
dependence, one can use the same integrated asymmetry, A' ud , as in Eq.flTBD except in this 
case the sums J^ 2 + and f£, 2 + J^/ 2 are automatically taken care of by the definition 
of 0o- Therefore Eq.(^) is still valid even in the case of l~l + qq or qqq'q' decay modes. Note 
that for these modes the asymmetries corresponding to Eq.(|T7D are not doubly suppressed 
by the small c v as in the purely leptonic modes. Unfortunately, A u ,± is not observable due 
to the lack of hadronic charge identification. However, one should keep in mind that if one is 
willing to zoom into the dependence of the amplitudes on the jet energies (or, equivalently, 
on the angles 9, or 9') then it may be possible to use energy identification to replace charge 
identification to construct an observable similar to A u ,d.- 

Note that for the CP conserving decay, the angular correlation we just described can 
also be used to detect the important cos 20 distribution even for the semi-hadronic modes, 
l~l + qq, or the hadronic modes, qqq'q', of the H — > ZZ decay. 

VI. H -» W+W- 

Once one understands how the CP-odd correlation can be decoded in the hadronic modes 
of H — > ZZ decay, it is easy to apply it to the process H — > W + W~ also. In this case one can 
write down a similar vertex with the Z fields replaced by the W fields. Since the coupling 
of the W boson to a fermion pair is always left handed, we just set cr = c' R = and 
cl — d L — 1. All substitutions are straightforward. We can simply add the superscript W 
to B, D, n, An, A u .d. and A' u d in Section IV to label the WW mode in the above formulas. 

In order to determine the decay planes of both W^ bosons, we cannot use the purely 
leptonic mode Ivvl, where too much kinematic information is carried away by the missing 
v and v. However, we can use the mixed modes Ivqq or Ivqq, where kinematics can be 
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fully reconstructed provided that we have good resolution on jets. Without identifying the 
nature of the leading quarks in jets, we are still able to define uniquely the 4>o as defined in 
Section V. This is good enough to measure the sin 20 dependence, which is related to A'^ d , 
a quantity generalized from the definition in Eq.(|TJ). In Fig. 3, we show A^ d , A'^ d versus 
Mh per unit of D w / B w . Note that, unlike the purely leptonic modes of — > ZZ case, 
A^ d is much lareger than A'^ d . Unfortunately, A^ d is not observable due to lack of charge 
identification. However, as commented earlier for the ZZ case, it may be possible to replace 
charge identification by energy identification to construct an observable similar to A^ d . 

For the purely hadronic mode, there is tremendous background from non-resonant con- 
tributions in the hadronic environment. The e~e + collider may be a cleaner machine in 
which the Higgs boson may be produced via the process e~e + — > ZH. Then, one can use 
the decay Z — > l~l + to tag the recoiling Higgs boson H which turns into a W + W~ pair and 
finally becomes 4 jets, whose CP character can be decoded using the angular correlation 
defined earlier. 

In addition, our study of the angular correlation can be very useful in investigating the 
general CP-odd or CP-even correlations in other reactions, for example, the hadronic final 
states of the process e + + e~ — > W + + W~ WM . 



VII. 

H -f it 

The possibility of detecting CP violation in H — > ti is very interesting in the sense that 
the corresponding CP violating couplings can occur at the tree-level. Let us look at the 
phenomenological form of the Yukawa interaction, 

C Yuka W a = _( mt / v )i( A p L + A*P R )W . (20) 

Here v = (y/2Gp)^^ ~ 246 GeV. The complex coefficient A is a combination of model- 
dependent mixing angles. The CP violating effect is proportional |T2|Jl3|] to Im A 2 = 



21m ARe A. In Higgs boson production at a high energy collider, one expects the CPT- 
even, CP violating effect to be observable in the angular asymmetry as before. It can in 
principle be tree-level effect and therefore can be very significant. However, in reality, its 
signal is harder to decode. The ti pair decay into W + bW~~b. If both W bosons decay lepton- 
ically, one faces the problem of identifying the Higgs boson event with the missing neutrinos. 
If a W boson decays hadronically one has to deal with its multi-jet final state. Assuming 
that one can identify the jets associated the W's then one can use the momenta of W's 
and b, b to define the decay plane of ti. Then the C-even angular correlation discussed in 
previous section can be used to decode CP violation. In hadronic colliders this is probably 
too hard to achieve. In a leptonic machine this may be possible only if the Higgs is produced 
in the futuristic NLC (such as EE500) or 77 colliders |15| . 



In this section, we will give the full differential form of the decay probability, which is 
useful for future experimental simulation. Without loss of generality, we look at the process 
H — > ti, with t —>■ bev and i —>■ bev. Formulas for other processes can be obtained by simple 
substitutions. Our result is based on the standard V— A coupling for the t decay . 
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]T \M\ 2 = 8g 8 m 2 t b ■ vb ■ VU2e ■ He ■ H - e ■ eH 2 )A] 

spin K 

+ [2e • (t - t)e ■ (t - t) - e ■ e(t - i) 2 ]A 2 R 

+ At^e^A^j^imt/vf^tViPw+Vw^ 2 . (21) 

For brevity, it is understood that we use the particle symbol to denote its corresponding 
momentum. The last four propagators, Vs, can be treated quite easily in the narrow- 
width approximation; simply replace the virtual momenta by their on-shell values, and the 
integrations over the squares of virtual momenta will give an overall factor 

f IX 71 \ 2 

\m t T t M W T W ) 

Inspecting Eq. (|2T|) , we know that e(e,e,t,t) is an appropriate CP-odd observable. Since 

e(e, e,t,i) = -\M H p e x p- e ■ (p t - ft) , 

in the H rest frame, CP information resides in the relative azimuthal angle between p e and 
p s with respect to the axis pt — pi- In contrast to Eq. (|TJ]) there is only one CP-odd angular 
distribution in Eq. (pT|) . This is because the top quark is a spin-| object and has only two 
helicity states. The corresponding CP-odd variable is thus proportional only to p+ — 
where ± represents the sign of the top helicities. 

In the future it may be possible to detect other similar CP-odd correlations in process 
such as 7 + 7— > t + i + H . 



VIII. 



e+e~ -> Z + H 

The analysis given above can be easily applied to the equally important process of the 
Higgs production as well. Here we shall simply investigate a typical one which is promising. 
The Higgs boson may be produced ]16[ through the Z bremsstrahlung in future high energy 



e~e + colliders and the CP information can be carried over by the lepton pair from the Z 
decay. To visualize the process of e~(pJ)e + (p + ) — > Z(P)H where Z — > Z(/c_)Z(/c + ), we sketch 
the reaction configuration in Fig. 4. The final on-shell Z boson is produced at a scattering 
angle 6. Its subsequent decay into a lepton pair is described by the polar angle 9 and the 
azimuthal angle (f) of I in the Z rest frame with its z axis opposite to P = k_ + k + . We focus 
our study in the case that the Z boson, produced simultaneously with the Higgs boson, is 
on-shell and we will use the narrow-width approximation to handle its decay process. 

The virtual gauge boson in the s-channel can be 7* or Z*. It has the momentum 
P* = p_ + p + . First, we write down the relevant form factors 

iM z * { p* tV ^ z{Pjri)H = tg H [Br] ■ 77* + D z *M H 2 e( V , 77*, P, -P*)] , (22) 

and 

iM r{P *^ z{Pin)H =ig H {D^Mv 2 e(r ] ,r ] *,P,-P*)] . (23) 
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The form factor D z * is related to D in Section II when one of the Z's becomes virtual. We 
put in a negative sign in front of P* so that the momentum flow is consistent with that 
defined in Eq.(2). The CP-odd term of D 1 * provides another source of CP violation. Other 
higher-dimensional CP-even terms are omitted in our discussion. The overall transition 
probability is 



4e 2 g H 



spin 



s-Ml 



1 



P 2 -Ml + iT z M z 



x lB< 



[ c t + C %)P- ■ k ~P+ ■ k + + 2c l c rP- ■ k +P+ ■ k 



+ Be(p_,p + ,k_,k^ 



c l c l + c r c r)(P- -k-+p+- k+) 



~ (c L c R + c R c L )(p^ -k + +p + -k_ 



/Ml 



(24) 



where the D form factors has been lumped into the following modified couplings, 



-2 



-2 



c\D z * 
c R D z * 



c L (l - M 2 z /s)D^ , 
c R (l - M 2 z /s)D^ . 



The Z couplings for the electron or the muon are cl = 
cr = tan6V- The differential cross section is 

/£|M| 2 A5dP 2 



\ + sin 2 6 W )/ (sin 6 W cos 9 



it- 



da 



d cos Qd cos 



8192vr 4 s 



(25) 
and 

(26) 



The kinematic factor As = 2\ph\/^/s is defined in the e~e + CM frame. The integration 
over P 2 can be simplified in the narrow width approximation. One just replaces P 2 by 
its on-shell value Mf and substitutes the corresponding propagator squared by 7i/(M z T z ). 
Under CP transformation, the configuration transforms as follows, 



fcos O, cos 9, 



(C0S(7T — O), COs(7T — 9), 



(27) 



If the differential cross section does not respect the symmetry of this transformation, CP is 
violated. The presence of the e term in Eq.(^) will produce such CP violation. One can 
integrate Eq.(p4|) over the polar angle 9 of I, and obtain a <fi distribution, 



da 



a 



'aM z (c 2 L + c 2 R )\ (c 2 L + c 2 R )B 2 S + AS 



d cos Qd(fi 64s 



6r ; 



sin 2 9w cos 2 9 



w 



Ml/s} 



(28) 



The second factor involving the Z width is just the branching fraction BF(Z — > It). We 
purposely separate the CP-even part S and the CP-odd part AE in the above formula. 
The CP-even part is 



Asin 2 + - 



(4 + sin 2 O cos 20) 



+ 



3ttM z 



Ml M 



H 



) sin O cos < 



.2 \2 
-R) 



(4 + 4) 2 



(29) 
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The CP-odd part consists of two terms, 

AS = (a 2 sin 2 9 sin 20 + a\ sin sin 0) A 5 . (30) 
Both terms of sin and sin 20 are CP violating. Their coefficients are 

Note that A = (1 + M| /s — Mjy / s) 2 — 4Mf / s. We can define the integrated asymmetries as 
before, 



Al H d =-\ rfcosG / -/ # - — ^r-j— = — tt °7"; 2/ , , (32) 
uA - aJ-i \Jo Jn J dcosBdcp 4(A + 12Mf/s) v ; 




3tt/2 ,2tt \ ^ a 




d cos Qdcf) 



tt(A + 12M|/s 



(33) 



In Fig. 5, we plot and A' u d per unit of D z * / B or D 1 * / B versus y/s for the case M fl = 80 
GeV. Both the CP-even and the CP-odd contributions to the cross section decrease with 
increasing s when s is far above the threshold. However the CP-even part decreases slightly 
faster and therefore both A ZI j and A' u d are increasing functions of s. In addition, A ZI j 
increases faster than A' u d because of the factor y/s in a\. The process rate is proportional 
to s _1 for large s, but since every reaction, including the background, is proportional to s^ 1 , 
we conclude that higher energy is preferable to detect CP violation provided one can get 
high enough luminosity. 

An analysis similar to what we have done here can also be applied to processes, like 
e~e + — > e~e + + (7*, Z*) + (7*, Z*) — > e~e + H. Of course even in e~e + — > ZH — > llbb, one 
can also analyze CP violation in the Higgs decay by zooming into the angular correlations 
of decay products from the final b quark jets. 



IX. HIGGS MODEL 



In renormalizable models the CP-odd couplings for both the H — > ZZ and the H — > 
W + W~ modes can be induced only at the loop level because they are higher dimensional 
operators. As the e symbol in Eq.(Q) only occurs through a fermion loop in perturbative 
calculations, it is natural to use the top quark as the internal fermion for its potentially 
large Yukawa coupling, see Eq. (p0|) . 

In this section, we will show the CP violating form factors D and D w in processes 
H — > ZZ and H — > WW based on the CP non-conserving Yukawa interaction in Eq. (p0|). 
The one-loop diagram of interest is shown in Fig. 6. We use the method of dispersion 
relations to find the form factor D for the process H — ► ZZ. The virtual mass s of the 
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Higgs boson is analytically extended beyond the on-shell value M#. In Ref. JIJ, we have the 
absorptive part, 



Im D(s 
with P 2 = 1- 



3a /2 2 P t 



47T ' < 4 



1 + 



K 



sin 2 # ^ 2 



K{3 2 



±PtPz v " 3 — 4A 
4m 2 /s, (3 2 z = 1 — AM^/s, and the logarithmic factor 

1 + /5 2 -2A/? Z 



ii 



(34) 



K = lOE 



1+/3 2 +2/% 



(35) 



The dispersive part is obtained by the Kramers-Kronig relation. 



Re D(s) 



-P 



71 



00 Im D(s') 



ds' . 



(36) 



The symbol P denotes the principal value of the singular integral. Fig. 7 shows typical sizes 
of Re D for various cases. It is of order of 10 -2 to in general. We must keep in mind 
that we only show the Higgs boson contribution in the perturbative regime. In scenarios 
where the Higgs bosons interact strongly, the CP violating effect in the colliders could be 
much larger. 

Similarly, the form factor D w for the process H — > W + W~ can be obtained with the 
following absorptive amplitude in Ref. [[!]]. 



Im D w (s) 



- — G F m t lm A-^-L(Pt,p w ^ 
4vr /% 



(37) 



with fiw = 1 — AMyy/ s and the corresponding logarithmic expression, 



L(x,y) 



x 2 — y 2 
2xy 



log 



x-y 



x + y 



(38) 



Fig. 8 shows typical sizes of Re D w for various cases. 



X. CONCLUSION 

We have shown that the angular asymmetry of the prompt lepton events originating 
in Higgs decay are sensitive to CP violation. Useful CP-odd variables are constructed in 
probing the dispersive parts of the CP violating form factors. We have demonstrated that 
the charge identification of jets is not necessary in order to look for CP signals in H — > ZZ 
or H —>■ W + W~. The expected event distributions which we give in their full differential 
forms are useful for future experimental simulation. Most of our approach to the problem 
is model independent, and so it is suitable for experimental search for the CP-odd form 
factors. 

While we were finalizing this manuscript, we came across two recent preprints on the 
subject of CP violation in the Higgs decay |L7,18 . Both papers have some overlaps with 



our work. We have noticed that our conclusion regarding the effect of having the identical 
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particles (H — > — > ////) differs from Ref. [|17[]. Also, some sign differences between 
our Eqs. (|l4lJ15|) and the corresponding ones in [17] can be attributed to the difference in 
definitions of <j). 

This work is supported in parts by the U.S. Department of Energy. DC wishes to thank 
L. Wolfenstein for many fruitful discussions. 
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FIGURE CAPTIONS 



Fig.l A typical decay configuration, H — > ZZ — > l(k-)l(k+)l' {p-)l' ■ 
Fig.2 A u . d ., A' u± versus M H /M Z per unit of D/B. 
Fig.3 AZ d ., A'Zd ver sus Mh/Mz per unit of D w /B w . 

Fig. 4 The configuration of e~(p_)e + (jo + ) — > Z(P)H where Z — > l(k^)l(k + ). Note that is 
defined in the rest frame of Z and is only represented schematically here. 

Fig.5 A z u H d and A' Z u H d per unit of D^/S or D^/B versus for M H = 80 GeV. 

Fig. 6 The one-loop diagrams which give the CP violating form factors D via the intermediate 
t quark. 

Fig. 7 Typical sizes of Re D for various cases. 
Fig. 8 Typical sizes of Re D w for various cases. 
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